Correlation-driven chiral superconductivity and chiral spin order in doped kagome 

lattice 
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We study the electronic instabilities of the Hubbard model in the 1/6 hole-doped Kagome lattice 
using the variational cluster approach. The 1/6 hole doping is unique in the sense that the Fermi 
level is at the von Hove singularity and the Fermi surface has a perfect nesting. In this case, a density 
wave is usually realized. However, we demonstrate here that the chiral d x 2_ y 2 +id xy superconducting 
state is most favorable when a small Hubbard interaction U(U < 3.0£) is introduced, and a scalar 
chiral spin order is realized at large U(U > 5.0t). Between them, a spin-disordered insulating state 
is proposed. 

PACS numbers: 74.20.-z, 75.10.-b, 71.10.-w, 71.27. +a 



The kagome lattice has recently attracted considerable 
interest due to its higher degree of frustration. Several 
possible states have been proposed for the Heisenberg 
model in this lattice, including the U{1) algebraic spin 
liquid(SL) fl, the valance bond solid [2| , and the gapped 
SL Q. Recently, the numerical study shows that its 
ground state is a singlet-gapped SL with signatures of Z 2 
topological order On the other hand, the anomalous 
quantum Hall effect and the topological insulator Q 
have also been demonstrated to exist when the electron 
filling is near the Dirac point at 2/3 electron density. In 
view of the rapid developments in the investigation of 
these exotic quantum phases, the question arises if other 
correlation-driven exotic quantum orders will be realized 
when the system is doped away from half-filling. We note 
that the Fermi level is at the von Hove singularity and 
the Fermi surface (FS) has a perfect nesting at 1/6 hole 
doping. By using the variational cluster approach(VCA) 
to the Hubbard model in the kagome lattice, we show 
that a chiral d x 2_ y 2 + id xy superconductiving (SC) order 
and a non-coplanar chiral spin density wave(SDW) can 
be realized at the 1/6 hole doping. 

Chiral superconductivity and chiral magnetic order 
are two distinctive phases of matter. They break both 
the time-reversal symmetry and the parity symmetry. 
The nontrivial topology of them can result in a wealth 
of fascinating properties, such as spontaneous quantum 
Hall effect [7|, |8( , unusual magnetoelectric properties [9| , 
and a quantized boundary current in magnetic field [lOj. 
Experimentally, the spin-triplet p-wave chiral supercon- 
ductivity has been found in Sr2Ru04 ll|. And, the 
chiral spin order was proposed to describe the mag- 



netic ordering in Mn monolayers on Cu(lll) surfaces [12 1 
and the nuclear spin ground state of a two-dimensional 
solid 3 He Recently, the chiral spin order was also 

proposed in the doped triangular lattice [Ljjthe py- 
rochlore lattice and the doped graphene 16| based 
on the mean-field analysis. In particular, the renormal- 



graphene with a perfect FS nesting. Here, we show that 
both the d x 2_ y 2 + id xy superconducting order and a non- 
coplanar chiral SDW order can be realized in the 1/6 
hole-doped kagome lattice by tuning the on-site Coulomb 
interaction U. 

The Hubbard model in the kagome lattice is defined as 
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ization group calculations [17|, LL8( show that the chiral 



d x 2_ y 2 + id xy superconductivity is favored in the doped 



where c i(T (ci a ) creates (annihilates) an electron with spin 
a on site i and — c\ a Cic (••) denotes the nearest- 
neighbor(NN) bond. Each unit cell of the kagome lattice 
contains three sites (labeled by a, f3 and 7 in Fig. 1(a)). 
Except for a flat band at the top of the energy band, 
the two dispersive bands are the same as those in the 
honeycomb lattice (Fig. 1(b)). From the density of states 
shown in Fig. 1(d), one can see that there are three von 
Hove singularities at — 2t, and 2t. That at 2t originates 
from the flat band, the other two originate from the sad- 
dle points at M point (see Fig. 1(b)). The Fermi levels 
at — 2t and correspond to the 1/2 and 1/6 hole dop- 
ing. In two cases, the FS forms a hexagon and displays 
a perfect nesting(see Fig. 1(c)). As the correlation effect 
at 1/2 doping is much weaker than that at 1/6 doping, 
we focus our study on the 1/6 hole doping here. 

VCA is a cluster approximation of the self-energy func- 
tional approach (l9| . which uses the rigorous variational 
principle 5tt(T,)/5T, = for the thermodynamic grand 
potential 57 to determine the physical self-energy S. It 
has been successfully applied to the problem of compet- 
ing phases in many strong correlation systems (20l423| . 
In VCA, the lattice is tiled into identical clusters, which 
make up of a reference system with the same two-body in- 
teraction as the original system but a different one-body 
part (including the added Weiss fields to study the sym- 
metry broken phases). Then the exact Green function 
G (and the self-energy £') for each cluster is calculated 
by exact diagonalization method. So, the short-range 
static and dynamical (within each cluster) correlations 
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FIG. 1: (color online) (a) Structure of kagome lattice and the 
12-site cluster tiling (enclosed by doted lines) used in our VCA 
calculation, ai and a2 are the lattice unit vectors, (b) The 
tight-binding dispersion along high symmetric directions(as 
illustrated in (c)). The dashed line is the Fermi level corre- 
sponding to the 1/6 hole doping, (c) The Billouin zone and FS 
for the 1/6 hole doping, bi and b2 are the reciprocal- lattice 
vectors. The dashed lines denote the FS and the vectors Q*' \ 
are the nesting vectors, (d) Density of state, (e) Weights of 
the contribution to FS from three inequivalent lattice sites a, 
P and 7 as represented by the colors. 



have been taken into account precisely. For any £ pa- 
rameterized as S'(t'), where t' represents the collection 
of the one-body terms, we have the grand potential [l9j : 

n[E'(t')] =0 / (t')+Trln[-G'(t / )] -Trln[-G(t')], (2) 

where O'(t') is the grand potential of the reference sys- 
tem and G(t') is the approximate Green function of the 
original system calculated through the cluster perturba- 
tion theory [24|]. Eq.@ is no longer a functional but an 
ordinary function of the variational parameters t', and 
the task of VCA is to find a stationary point of this func- 
tion, dVt(t')/dt' = 0. In our calculation, the 12-site clus- 
ters^ enclosed by the dotted lines in Fig. 1(a)) and the 
open boundary conditions for the clusters are used. 

Before presenting our numerical results, let us first dis- 
cuss the possible SC symmetries at 1/6 hole doping. We 
note that at each saddle point M the site-contribution 
to FS comes only from one of the three inequivalent lat- 
tice sties, as shown in Figl.(e). Considering the effect 
of the van Hove singularity at Af-point, it is expected 
that the favorable Cooper pairings will be made from 
two electrons belonging to the same sublattice. For the 
kagome lattice, there are six coordinates in the same sub- 
lattice for each site (Fig. 1(a)), and it gives rise to three 
different bonds. Considering that the paring intensi- 
ties are the same along all bonds, and the phase dif- 
ferences 9 and </> with respect to that along the a2 di- 
rection exists for the other two directions, we will get 
the following pairing functions: (i) A(k) = cos(k x ) + 
e w cos(fc x /2 + V3k y /2) + cos(fc 2; /2 - V3k y /2) for the 
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FIG. 2: (color online) (a)-(d) Phase(±) of the SC gap func- 
tions for the four possible pairing symmetries. The dotted 
lines denote the gap nodes, (e) Scaled grand potential Q as a 
function of Weiss field A for different U in the d x 2_ y 2 + id xy 
channel, (f) Scaled Q as a function of A for different pairing 
symmetries at U — t. 



V3ky/2) + s\x\{k x /2 - V?>k y /2) for the spin-triplet 
pairing. The most natural choices for the phases are 
{6,<t>) = (0,0) and (2tt/3, -2tt/3). So, we have: (i) 



the d x 2_ y 2 + id xy wave for {6,<f)) — (27r/3, — 27r/3), 
A d+td (k) = A [cos(fc x ) - cos{k x /2)cos{^/Zk x /2) + 
sm(k x / 2) sin(v / 3fc a; /2)]; (ii) the extended s wave 
for (0,0) = (0,0), A s (k) = A [cos(fe x ) + 
2 008(^/2) cos(v / 3fcx/2)]; (iii) the p x + ip y wave for 
{9,cp) = (2tt/3, -2tt/3), AP+^(k) = & [sm{k x ) - 
sin(fcc/2) cos(\/3fca;/2) + j-\/3 cos(fc :c /2) sin(v / 3fc 2; /2)]; and 
(iv) the / wave for (0, 4>) = (0, 0), A'(k) = A [sin(fc a; ) + 
2 sin(fc x /2) cos(v / 3fc a; /2)]. The phases and the gap nodes 
for the four pairing symmetries are shown in Fig.2(a)-(d). 



With the above pairing functions, we now introduce 
the Weiss fields H' = AY!/,.,h e i9 «A, s f + h.c, where 
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spin-singlet pairing, (ii) A(k) = sin(fc 2 
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sia(A; x /2 + 



c it c jl T c il c jt f° r the singlet/triplet pairing, to 
test the possible pairing orders. Here, ((••)) denotes the 
bonds linking the sites in the same sublattice and Oij 
is given according to the rule discussed above. Fig. 2(e) 
presents the results of fl — Oo as a function of A for 
the d x 2_ y 2 + id xy pairing symmetry at various U. f2o 
is the grand potential in the zero Weiss field. For 
OAt < U < 3i, we find that there is a minimum at fi- 
nite A, at the same time it satisfies <9f2(A)/<9A = 0. For 
U < OAt, a monotonic increase with A occurs. While, 
though there is also a minimum at finite A for U > 3i, 



3 



we find that the derivative at this minimum does not ex- 
ist, because Q is not smooth at the point. Thus, it is 
not a stationary point of the self-energy functional, ac- 
cording to the variational principle of VCA [l9j]. So, the 
d x 2_ y 2 + id xy pairing order is realized only in the region 
OAt < U < 3t. We have also checked the stationarity 
of other three SC pairing orders. The typical results are 
shown in Fig. 2(f) for U = t. Except for the d x 2_ y 2 + id xy 
pairing, all other three pairings exhibit a monotonic de- 
crease with the Weiss field. Therefore, they are not stable 
solutions. 

In the system with a perfect nesting FS, one will ex- 
pect naturally an instability to the density wave which 
may compete with the SC order. In the presence of 
the on-site Hubbard U, the most likely density wave 
would be a SDW. A well known example is that the 
staggered SDW order occurs for the half-filled Hubbard 
model with a perfect FS nesting in the square lattice, 
where the nesting wave vector (±7r, ±7r) coincides ex- 
actly the real-space translation symmetry of two unit 
distances for the staggered SDW. Considering the geom- 
etry of the kagome lattice and consequently the nesting 
property for the 1/6 hole doping as shown in Fig. 1(c), 
the staggered SDW will not be favored here. Then, 
what is the kagome-lattice counterpart of the staggered 
SDW occurring on the half-filled square lattice? At 
1/6 hole doping, the nesting FS has the nesting vec- 
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tors Q^ c (or Q„ b c which connects to 
reciprocal-lattice vectors) as indicated in Fig. 1(c). With 
these nesting vectors, two distinct magnetic orders are 
expected. One is the non-coplanar chiral spin order asso- 
ciated with b c as shown in Fig. 1(c), which is sim- 



FIG. 3: (color online) (a) Non-coplanar arrangement of 
spins(scalar spin chirality). Si-h indicate the four direc- 
tions of local magnetic moments, (b) Coplanar arrangement 
of spins (vector spin chirality). The dashed lines with arrow in 
(a) and (b) are the magnetic translation vectors, (c) Scaled 
grand potential Q as a function of the Weiss field h for the 
spin order shown in (a) at various U. The local minima are 
indicated by arrows, (d) Scaled Q as a function of h for the 
spin order shown in (b). The inset of (d) shows the energy dif- 
ferences between the two spin ordered states and the normal 
phase with a zero Wiess field, respectively. 



To test the favorable SDW, we introduce the following 
Weiss field: H' = hJ2i V e «j ' Smj, where is the orien- 

J 



ilar to that proposed in the doped triangle [14{ and 
honeycomb [16( lattices, where the four local spins di- 
rect along the normals to the faces of a regular tetra- 
hedron(Fig.3). Thus, the spin orientations are: ej a = 
l/\fZ[e x cos(Q* ■ Rj) +e y cos(Q b • Rj) +e z cos(Qj ■ Rj)], 



tation of magnetic moment and Si v = J2a<r' c iria' r ^' °m"< 
with r the Pauli matrixes, i the unit cell index and n the 
site index. Fig. 3(c) and (d) show the results for ft — Qq 
as a function of the Weiss field h for the scalar chiral 
spin order and the vector chiral spin order, respectively. 
We find that a local minimum exists only for U > 5.5t 
in both cases. Therefore, in the region OAt < U < 3i, 



1/V3[e x cos(Q* - Rj) — e y cos(Q b ■ Rj) — e z cos(Qj • no SDW state will compete with the chiral d x 2_ y 2 + id x 



e»/3 _ 

Ri)] and e. n = l/y/3[-e x cos(Q* ■ Rj) - e y cos(Q£ ■ Rj) + 
e 2 cos(Qc - Rj)], where Rj denotes the position of the unit 
cell i. In this state, the resulting scalar spin chirality, 
(K s ) = (S Q • (S/3 x S 7 )) 7^ in each triangular plaquette, 
breaks both the time-reversal and the parity symmetries. 
The other is the coplanar vector chiral spin order asso- 
ciated with 0,1 bc , as shown in Fig. 3(b). In this case, 
the spins in each triangular plaquette orient at 120° to 
each other. The vector chirality for each triangle can be 
defined as K„ = (2/3-s/3)(S a x Sp + Sp x S 7 + S 7 x S a ), 
which is parallel to the z axis with amplitude +1 or —1. 
As shown in Fig. 3(b), the vector chirality arranges as a 
"stripe" phase, in which each chirality arranges in a stripe 
and the two stripes are staggered with each other. The 
scalar chiral spin order and the vector chiral spin order 
have different translational symmetries in real space as 
indicated in Fig. 3(a) and (b), where the magnetic trans- 
lation vectors are shown as the dotted lines. 



pairing state. Fig. 3(c) and (d) show that both spin or- 
dering states are more favorable compared to the spin 
disordered state. So, we need to compare their energy 
gain which is given by the energy difference between the 
ordered solution and the normal solution found by sup- 
pressing the Weiss fields. The energy density E = CI + fin 
as a function of electron density n can be obtained for 
a solution with the functional ft. The calculation of the 
energy gain is performed for several values of the chemi- 
cal potential fi until the density n is close enough to 1 /6 
doping, and the cluster's chemical potential fi' is also 
used as a variational parameter to guarantee the ther- 
modynamic consistency 21] . The inset of Fig. 3(d) shows 
the energy gain for the two spin ordered solutions as a 
function of U. It shows that the scalar chiral spin order 
state is more favorable than the vector chiral spin order 
state for U > 5.5t. With decreasing U, the difference of 
the energy gain for the two spin-ordered states decreases. 
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FIG. 4: (color online) (a) Scaled grand potential Q as a func- 
tion of Weiss field A for the conventional CDW order (see 
text) at various U. (b) Spectral functions for U = At along 
the high symmetric directions(see Fig. 1(c)). (c) Qualitative 
phase diagram of the Hubbard model in the kagome lattice 
for 1/6 hole doping. 



Near the threshold U for the occurrence of the spin order 
states, the two states are nearly degenerate. 

To check the possible existence of the charge-density 
wave(CDW) in this system, H' = A -(/ Q n iQ cos(Q^ • 
r») + fpriip cos(Q^ • + / 7 n i7 cos(Qj • r,)) is used, with 
i the unit cell index and (f a , fp, / 7 ) = (1,0,-1). The 
dependence of il — firj on the Weiss field A is presented in 
Fig4.(a). We find that SI — ilo decreases monotonously 
with A for various U from U = 2 ~ Qt. Therefore, the 
CDW does not exist in this region. Thus, in the region 
OAt < U < 3t, the only stable state we found is the 
chiral d x 2_ y 2 + id xy pairing state, while for U > 5.5t 
the most favorable state is the scalar chiral spin order 
state. In the region between them, we do not find ev- 
idence of the SC state, spin order state and the CDW 
state. From the spectral function of single-particle exci- 
tations, we find that a gap occurs at the Fermi level as 
presented in Fig. 4(b) for U = At (a non-vanishing spec- 
tral weight can be seen in the gapped region because a 
finite Lorentzian broadening is used in the numerical cal- 
culations). So, the system is in fact an insulator. We 
would suggest that the system is a magnetic disordered 
insulator (MDI). The possible candidates may be the spin 
liquid, the spin glass or the valence-bond solid, but the 
method used here can not identify or distinguish them. 

Finally, we summarize the results obtained above in 
the phase diagram shown in Fig4.(c). We note that the 
spin- 1/2 kagome lattice has been realized in Herbert- 
smithite ZnCu3(OH) 6 Cl2 [1^, (2(| and its isostructural 



Mg-based paracatamite MgCu 3 (OH) 6 Cl 2 [271. Also > tne 
kagome lattice has been simulated experimentally in 
ultra-cold atoms [28| ■ So, we suggest that the theoretical 
predictions presented here may be probed after doping 
these compounds or by implementing an optical lattice 
in ultra-cold atoms where the Hubbard interaction U can 



be tuned continuously. 

In summary, we have studied the Hubbard model in 
the 1/6 hole doped kagome lattice using the variational 
cluster approach. We find that the chiral d x 2 -y 2 + id xy 
superconducting state is most favorable at a small Hub- 
bard interaction U(0.At < U < 3t), and the scalar chiral 
spin order is realized at the large U(U > 5.5t). Between 
them, a magnetic disordered insulating state is proposed. 
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